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ABSTRACT 



We present the result of investigations into two theories to explain the star for- 
mation rate-density relationship. For regions of high galaxy density, either there are 
fewer star forming galaxies, or galaxies capable of forming stars are present but some 
physical process is suppressing their star formation. We use HI Parkes All Sky Survey's 
(HIPASS) HI detected galaxies and infrared and radio fluxes to investigate star forma- 
tion rates and efficiencies with respect to local surface density. For nearby (vel< 10000 
km s _1 ) HI galaxies we find a strong correlation between HI mass and star formation 
rate. The number of HI galaxies decreases with increasing local surface density. For HI 
galaxies (1000<vel<6000 km s _1 ) there is no significant change in the star formation 
rate or the efficiency of star formation with respect to local surface density. We con- 
clude the SFR-density relation is due to a decrease in the number of HI star forming 
galaxies in regions of high galaxy density and not to the suppression of star formation. 

Key words: methods: galaxies:formation - evolution - spiral - radio lines: galaxies 
- radio continuum: galaxies - infrared: galaxies. 



1 INTRODUCTION 



iDressler (1980ft . finds a well-defined relationship between 
local galaxy density and morphology for galaxies in 
clusters, the 'morphology-density re la tion'. Several late r 
studies, e.g. lHashimoto et al. C 1998ft. iGoto et al. (2003ft . 
ISodre fc Mateus (2003ft and lNuiiten et a! (2005ft use large 
optical surveys and 2 or 3D local galaxy density to in- 
vestigate the Morphology-density relation. They confirm 
IDressler (1980ft 's cluster results: as local galaxy density in- 
creases, there is a decrease in the relative number of spiral 
galaxies and an increase in the number of elliptical and SO 
galaxies. 

Two studies investigating the morphology-density rela- 
tion with respect to the star formation rate (SFR) and lo- 
cal galaxy density (her eafter the S FR-density relation) are 
iLewis et al. (2002ft and lGomez et al. (2003ft . These studies 
are based on optical observations with galaxy SFRs calcu- 
late d from Ha fluxes. T o de termine the local galaxy den- 
sity |G6mezeiaO20^ use the lo- 
cal surface density. ILewis eta!^^002F) uses the SFR for a 
given density range for clusters compared to field galaxies. 
iGomez et al. (20 031 uses two samples, one sample is sepa- 
rated into field and group/clusters gal axies and the other 
is the Sloan Digital Sky Survey. Both ILewis et al. (2002ft 



and lG6mez et aL (2003) conclude the SFR for a given den- 
sity range decrease s as the local surface density increases. 
iBalogh et al. (2004ft investigate the relationship between the 
distribution of galaxies and local density and find that star- 
forming and quiescent (elliptical) galaxies vary with local 
density. However, they find the 'distribution of Ha equiva- 
lent widths' of the star-forming population does not depend 
strongly on environment. 

Our work is based on the HI Parkes All Sky Survey Op- 
tical Catalogue (HOPCAT), dPovle et al. 2005ft . HOPCAT 
is, at the present time, the largest optically matched HI ra- 
dio detected catalogue of the whole southern sky. For the 
first time we investigate the SFR-density relation using a 
large sample of galaxies containing the basic fuel for star 
formation, HI. For regions of high galaxy density we ask, 'is 
the decrease in SFR due to fewer star forming galaxies or are 
galaxies capable of forming stars present but some physical 
process is suppressing their star formation?' 

When investigating the SFR-density relation, vari- 
ous 3D methods to determine t h e local galaxy den- 
sity are used iAllam et al. 19991 : iMartmez et al. 20021 : 
ISodre fc Mateus 2003ft ' which require redshifts. At the time 
of writing, no redshift catalogue of the whole of the 
southern sky exists, hence we use the local surface den- 
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sity, Sip, similar tolDressler ( 1980ft . lLewis et al. (2002) and 
iGomez et al. (2003ft . 

To calculate the SFR. iHookins et al. (200lft . 
Cardiel et al. ( 2003ft . iLewis et al. (2002ft and 

Gomez et al. (20031) use Oil, OIII, Ha, H/3, H<5 and 
UV detected fluxes, however these bands require correction 
for dust extinction. Infrared and radio detections, both 
indic ators of current star formation, are no t affected by 
dust (iHopkins et al. 20031: ICardiel et al. 2003h and studies 
show a correlation exists between IR a n d 1.4-G Hz SFRs 
jHopkins et al. 20021 : ICardiel et al. 20031 lYun et al. 200 it) . 
In our study we use 843-MHz and 1.4-GHz radio, and the 
12, 25, 60 and 100 (jm infrared (IR) detected fluxes to 
calculate the SFR per HI galaxy. 

The star formation efficiency (SFE) is used to determine 
the efficiency with which the HI present is converted into 
stars per year. From the HI galaxy SFR, SFE and Eio we 
explain the reason behind the SFR-density relation. 

In Section |2] we describe the process taken to deter- 
mine the local surface density for each HOPCAT target 
galaxy. Section [3] details the resources we use to obtain the 
fluxes and the equations to calculate the SFRs and SFEs. 
We present our results in Section[I]and discuss these results 
in Section^] Section|S|summarises our work. We use H =71 
km s _1 Mpc -1 throughout. 



2 LOCAL SURFACE DENSITY 

To investigate the SFR-density relation we require a mea- 
sure of local galaxy density. Using the positions for optically 
matched HOPCAT galaxies and two optical background cat- 
alogues, the local surface density for each galaxy is cal- 
culated. In the following sections we introduce the target 
galaxies from the largest optically matched HI survey of the 
whole southern sky, HOPCAT. We also discuss the back- 
ground catalogues, the S uperCOSMOS Galaxy Catalogue 1 
Hamblv et al. 2001allblcft and the 6dFGS input catalogue 
Jones et al. 2004ft that are used, along with HOPCAT, to 
calculate the local surface density, Eio. 

2.1 HOPCAT Target Galaxies 

Our base sample is optically matched HOPCAT galaxie s 
from match categories 60, 50, and 40, iDovle et al. (2005ft . 
These Hi-optical galaxy match categories are; 60s - matched 
using independent velocity confirmation, 40s - matched us- 
ing independent velocity confirmation but multiple galaxies 
with similar velocities are within the 15x15 arcmin image 
and 50s - matched without independent velocity confirma- 
tion but the chosen galaxy is the only galaxy visible within 
the image (for a more detailed desc ription of HOPCAT and 
the match categories please refer to lDovle et al. (2005ft ). As 
these are the most reliable match categories we use their 
velocities in the local surface density calculations. In the 
subset of HOPCAT analysed in this paper the match cat- 
egory breakdown is 72 per cent 60s, 24 per cent 40s and 4 
per cent 50s. 

1 SuperCOSMOS Galaxy Catalogue is supplied directly from 
Mike Read at Royal Observatory Edinburgh 
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Figure 1. SuperCOSMOS Galaxy Catalogue sky coverage. Note 
the galactic plane (|6| < 20°), the Large and Small Magellanic 
Clouds and Dec>+3.25 are cut areas. Any HOPCAT target 
galaxy with less than 50 per cent of the 10th neighbour area 
overlaps a cut area, have corrected densities. Any HOPCAT tar- 
get galaxy within a cut area or with greater than 50 per cent of 
its 10th neighbour area within a cut area, are removed from our 
target lists. 

The basis of HOPCAT is the HI All Sky Surve y 
catalogue (HICAT), llMever et al. 20041: IZwaan et al. 2001) . 
This is a flux limited catalogue due to the sensitivity of the 
Parkes 21-cm multibeam receiver. Due to this, and the mag- 
nitude limits of the two background galaxy catalogues (see 
Section 12.21 , we only use galaxies between 1000 and 6000 
km s _1 for the SFR, SFE and Eio analysis. 

2.2 Background catalogues for Density 
Calculations 

We use two whole southern sky catalogues, SuperCOSMOS 
Galaxy Catalogue and The 6dFGS' 2MASS input catalogue 
to calculate the local surface density. 

The SuperCOSMOS Galaxy Catalogue (hereafter 
SCosGC), a subset of the SuperCOSMOS Sky Survey which 
contains only visually confirmed galaxies, is supplied di- 
rectly by Mike Read from the Royal Observatory Edinburgh. 
SCosGC contains approximately 107000 objects, has a faint 
magnitude limit of Bj < 17.0 mag and covers the whole of 
the southern sky up to +3.2° Dec. The areas around the 
galactic plane (|6| < 20°) and the Large and Small Magel- 
lanic Clouds are not included, Fig.0 Read (Private Commu- 
nication) compares SCosGC with the APM catalogue and 
finds a completeness > 95 per cent. The resulting densi- 
ties calculated using SCosGC are referred to as HOPCAT- 
SCosGC. 

The 6dFGS' 2MASS input catalogue (hereafter 
6dF2MASS) is o btained from the 6dFGS web site 
(http://www.mso.anu.edu.au/6dFGS/ 1. These objects are 
detected in the near-infrared with a K-band faint magni- 
tude limit of K < 12.75 mag. 6dF2MASS covers the whole 
of the southern sky up to -00° Dec. and contains approx- 
imately 121000 objects with the galactic plane (\b\ < 10°) 
not included, Fig. [5] A conversion of <Bj-K>=4 mag, cal- 
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Figure 2. 6dF2MASS Galaxy Catalogue sky coverage. Note the 
the galactic plane (|6| < 10°) and Dec>0 are cut areas. Any 
HOPCAT target galaxy with less than 50 per cent of the 10th 
neighbour area overlaps a cut area, have corrected densities. Any 
HOPCAT target galaxy within a cut area or with greater than 50 
per cent of its 10th neighbour area within a cut area, are removed 
from our target lists. 



culated from galaxies common to SCosGS and 6dF2MASS, 
is used to convert 6dF2MASS' K to Bj mags to standardise 
all magnitudes in this work. The densities calculated using 
6dF2MASS are referred to as HOPCAT-6dF2MASS. 



2.3 Density Calculation 

We use local surfa ce de nsity, Sip, similar tol5ressler 
iLewis et al. (20021) and lGomez et al. (2 003). to calculate the 
local galaxy density. The area of density is calculated from 
a Mpc radius (r) between the position of the HOPCAT tar- 
get galaxy and the 10th nearest neighbouring galaxy us- 
ing SCosGC and 6dF2MASS galaxy positions. As this area 
(hereafter 10th neighbour area) includes a total of 11 galax- 
ies, the Eio density is calculated as follows: 

E 10 (M P c- 2 ) = (1) 

As Eio is based on a 2D view of the sky, 
fore and back ground galaxie s can falsely inflat e the 
density value. iDressler (198Cf) . ILewis et al. (2003) and 
iGomez et aT72 003l use a faint absolute magnitude limit 
equivalent to Mb .=-17 mag to eliminate the effect from 
background galaxies. As we are using nearby HI detected 
galaxies, where fainter galaxies are detected, we require a 
fainter absolute magnitude limit as well as a foreground 
limit. We use Mb =-21 and M=-17 mag for arbitrary bright 
and faint absolute magnitude limits respectively. The cor- 
responding bright and faint apparent magnitude limits are 
calculated using the HOPCAT target galaxy's velocity and 
only galaxies within this range in SCosGC and 6dF2MASS 
are used to determine the 10th nearest neighbouring galaxy. 

Edge-affected regions occur near the areas cut from 
SCosGC and 6dF2MASS, Figs.[]]and|2| To reduce this prob- 
lem we use a number of culling steps. Any HOPCAT target 
galaxies in edge-affected regions, where less than 50 per cent 



of the 10th neighbour area overlaps a cut area, have cor- 
rected densities. Any HOPCAT target galaxy within a cut 
area or with greater than 50 per cent of its 10th neighbour 
area within a cut area, is removed from our target lists. Also 
any galaxies, mainly those with small velocities, having fewer 
than 10 neighbouring galaxies are removed from the tar- 
get lists. HOPCAT-SCosGC and HOCAT-6dF2MASS tar- 
get lists contain 2126 and 2603 galaxies respectively for 
1000<vel<6000 km s" 1 . 



3 STAR FORMATION RATES & EFFICIENCY 

To investigate star formation and hence the SFR-density 
relation, the HI galaxies SFRs and SFEs are calculated. The 
SFR is the rate of star formation in solar masses per year. 
SFE is the efficiency at which the HI present in galaxies is 
converted into stars per year and is independent of distance. 
We use the two HOPCAT target lists, HOPCAT-SCosGC 
and HOPCAT-6dF2MASS, to search for IR and radio fluxes 
to calculate the SFR and SFEs. Hereafter the 4 data sets will 
be referred to as shown in Tabled 

In the following section we discuss the IR and radio data 
resources to obtain fluxes and the equations to calculate the 
SFR and SFE. We also discuss the difference between the 
IR and radio SFRs. 



3.1 SUMSS and SN6I Catalogues 

To calculate SFR and SFE we utilise 2 sources of radio 
data, The Sydney U niversity Molonglo Sky Survey (SUMSS) 
iMauch et al. (20031) and a catalogue of SUMSS, NRAO 
VLA Sky Survey (NVSS), 6dFGS and IRAS da ta (hereafter 
SN6I) compiled by Tom Mauch jMauch 20051) . SUMSS is 
detected at 843-MHz and has a sky coverage of -90<Dec.<- 
30. NVSS is detected at 1.4-GHz and has a sky coverage of 
-40<Dec.<0. 

To determine the correct SUMSS and SN6I object corre- 
sponding to the HOPCAT target galaxies, we use the follow- 
ing analysis to find the optimum search radius for match- 
ing the respective catalogues. We compare the number of 
SUMSS and SN6I galaxies matched to HOPCAT using the 
original galaxy positions with those matching a 'random' 
catalogue made with the original positions offset by 0.25 
degs. The optimum search radius we use is where the num- 
ber of (original position) matches converges to the 'random' 
level denned by the offset positions. The optimum search ra- 
dius is 30 arcsec with a reliability of a match and complete- 
ness (the percentage of possible genuine matches missed) of 
80 and 1.24 per cent respectively. Final target galaxy num- 
bers with SUMSS and SN6I fluxes are listed in Table H 

The equations we use to calculate the SFR jCram 19981) 
and SFE, from the SUMSS and SN6I fluxes, are: 



SFRnajioiMQyr- 1 ) = 



LlA-GHz 



2.3 x 10 22 WHz 



SFE Radio (yr 



SFRiA-GHz 



(2) 



(3) 



The uncertainty for the 1.4-GHz radio SFR equation 
is a factor of 10 ' 2 , as discussed bv ICram (19981) . and is 
associated with the luminosity calculation. There are also 
random uncertainties in the SUMSS flux as listed in the 
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SUMSS catalogue, and an uncertainty of 0.1 per cent in 
SN6I fluxes as described in lMauch (20051) . The total ASFR 
and ASFE are: 



ASFR = SFR 



ASFE = SFE 



AF 2 
F 



2AVel 2 ASFR eq 2 



Vel 



SFR 



AF 2 
F 



AS„ 

Sini 



; 2 ASFR 2 
+ SFR ' 



(4) 



(•>) 



where F is the flux density, Si„t is the HI integrated flux and 
ASFR eq is the uncertainty in the 1.4-GHz SFR equation. As 
the SFR equation is for 1.4-GHz detected fluxes, the SUMSS 
fluxes are adjusted usi ng a mean rad io spectral index cal- 
culated as described in lMauch (20051) . consistent with pre- 
vious spectral index calculations for frequ e ncies below 1.4- 
GHz jMauch et al. 2o5H lOort et al. 19881 iHunstead 199lt 
iDe Breuck et al. 2000l) ." The radio spectral index is defined 
as a power- law: 



F 1A -GHz = F 8 43-MHz (1400/843) c 



(6) 



where F is the flux density, a— -0.7 is the radio spectral index 
and (1400/843) is the frequency ratio. 



3.2 IRAS 

We obtain our second source of fluxes from the Infrared 
Astronomical Satellite (IRAS) 2 . We use the IRAS Galax- 
ies and Quasars Catalogue (IRASGQCat) which covers the 
whole sky. The same method we use to determine the search 
radius for the SUMSS and SN6I objects is applied. The op- 
timum search radius of 60 arcsec between the HOPCAT and 
IRAS objects with a reliability and completeness of a match 
of 97 and 2.7 per cent respectively. Using this radius we 
search the IRAS web based batch search facility, GATOR 
{ http: / /irsa. ipac.caltech.edu /cgi-bin/Gator/nph-dd ) . using 
our HOPCAT target galaxy positions to obtain flux mea- 
surements. Final target galaxy numbers with IRAS fluxes 
arc listed in Table H 

The IRAS 12, 25, 60 and 100 /im fluxes are combined to 
calcul ate a total luminos ity for each galaxy from which the 
SFR dCardiel et al. 20031) and SFE are calculated as follows: 

FTot.m(ergs cm 2 s) = 1.8 x 10 11 
(13.48F 12M + 5.16F 25M + 2.58F 60m + F 10 o M ), 

(7) 

SFR T ot.in( Moyr" 1 ) = 4.5 x lQ- iA L To t.iR (ergs s" 1 ), (8) 
SFRrot .IR 



SFETot.miyr- 1 ) 



Me 



0) 



We use only the uncertainty, ± 30 per cent, ass oci- 
ated with t he SFR equation from ICardiel et al. (20031) as 
discussed in lKennicutt (19 981. as any random uncertainties 
are insignificant in comparison. 



2 This research has made use of the NASA/ IPAC Infrared Sci- 
ence Archive, which is operated by the Jet Propulsion Laboratory, 
California Institute of Technology, under contract with the Na- 
tional Aeronautics and Space Administration. 



Table 2. Galaxy samples. 



Fluxes 


HOPCAT-SCosGC 


HOPCAT-6dF2MASS 


SFR vs M HI 


(vel < 10000 km s" 1 ) 




IRAS 


1405 


1405 


SUMSS 


340 


340 


SN6I 


221 


221 


Sio Plots 


( 1000<vel<6000 km s" 1 ) 




IRAS 


923 


1180 


SUMSS 


172 


275 


SN6I 


167 


203 



3.3 Comparing IR and 1.4-GHz Radio based 
SFRs 

The SFRs from IR and 1.4-GHz radio are based on two dif- 
ferent flux producing mechanisms. The IR flux is a measure 
of the dust and gas heating during the star formation process 
and 1.4-GHz radio flux is a measure of sychrontron emission. 
Both are measures of current star formation and although 
they are from different mechanisms, the dust transparent 
nature of both measurements and the linear correlation be- 
tween 1.4-GHz radio and IR-based SFR, ensure both are 
equally reliable measures of current star formation. However 
as the uncertainties for IR are ± 30 per cent, the 1.4-GHz 
SFRs are more reliable. 



4 RESULTS 

In this section we present the results for our HI mass (Mhi) 
and SFRs investigations for HOPCAT galaxies. We also dis- 
cuss the difference between HI radio and optically detected 
galaxy samples and how this relates to SFRs. The SFRs and 
SFEs with respect to £io results are also presented. Please 
note: all scatter plots use least-squares line of best fit, cal- 
culated using the logs of Mhi, SFR or SFE, with ±3<r slope 
uncertainty lines. 

4.1 HI mass-SFR Relationship 

As HOPCAT is the largest optically matched HI radio de- 
tected catalogue, we comprehensively investigate the rela- 
tionship between Mhi and SFR. 

We use the integrated flux (Sint) and the velocity (vel) 
to calculate the Mhi as follows: 



Mhi(M q ) = 2.356 x 10 5 



(vel\ 2 
\ ti J 



(10) 



The AMhi, resulting from the velocity and S, 



l|Zwaan et al. 20041) . is 

Avel = 1.0 x lO 4 ^ 2 (mJy) + 5, 



ASim — 0.5Sl n t, 



1 2 Avel 2 AS lnt 2 
AM hi = M H i\l — + 



vel 



Sint 



(11) 

(12) 
(13) 



The Avel d epends only on the peak flux S p 
JZwaan et al. 20041) . 
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Table X. Data Set Names. 



Fluxes HOPCAT-SCosGC HOPCAT-6dF2MASS 



IRAS HOPCAT-SCosGS-IR HOPCAT-6dF2MASS-IR 

SUMSS & SN6I HOPCAT-SCosGS-1.4GHz HOPCAT-6dF2MASS-1.4GHz 



We calculate Mhi for all HOPCAT target galaxies with 
SFRs out to 10000 km s" 1 (IR: 1405 galaxies and 1.4-GHz: 
561 galaxies). Fig. [3] shows the relationship between Mhi 
and SFR with least-squares lines of best fit, as described 
earlier. AMhi are shown in Fig. |3] as; IR, red (min) and 
black (max); 1.4-GHz, Sumss, red (min) and black (max) 
and NVSS, magenta (min) and cyan (max). It is seen for 
both IR and 1.4-GHz, the Mhi increases with increasing 
SFR. 

However the flux-limited nature of our data is seen in 
Fig. [3]in the colour/symbol velocity groupings. For a given 
velocity only galaxies down to the sensitivity of the receiver 
are detected. To reduce this effect we analyse Mhi versus 
SFR for 1000 km s" 1 velocity ranges and fit each with a line 
of best fit as described above. These results also show an 
increase in Mhi with increasing SFR. 

Though Fig. |S|and the velocity separated analysis show 
Mhi increases with increasing SFR, the flux-limited nature 
of our sample may still affect results. To reduce the effect we 
use a volume limited subsample of our HOPCAT target list. 
We use only galaxies out to 4000 km s~ with M hi greater 
than 3 x 1O 9 M (see fig. 2 of lZwaan et al. f 20051) for velocity 
and Mhi limits). Fig.0]shows the Mhi-SFR relation for the 
volume limited sample. Again we use least-squares line of 
best fit. The Mhi still increases with increasing SFR. The 
gradients (m) in Fig. El are 0.6±0.02 and 0.4±0.02 for IR 
and 1.4-GHz respectively and m=0.3±0.02 for the volume 
limited sample, Fig. 01 

Though the gradient for the volume limited sample is 
less steep, again there is a clear increase in Mhi with increas- 
ing SFR. We conclude, there is a strong correlation between 
a galaxy's Mhi content and SFR. 



4.2 Comparing HI radio and Optical galaxy 

samples with respect to local surface density 

We use the 6dFGSs Second Data Release (hereafter 
6dFGS2DR) , obtained using their web access page 3 and our 
HOPCAT target galaxies to compare HI and optically de- 
tected galaxy samples for the southern sky. We have limited 
the 6dFGSDR2 data to contain only galaxies with the same 
velocity limit as our HOPCAT target list (1000<vel<6000 
km s _1 ). 

The left panel of Fig. |S] illustrates the difference be- 
tween HI and optically detected galaxies with respect to 
Eio. Though the 6dFGSDR2 sample (3901 galaxies) is larger 
than the HOPCAT target list (2602 galaxies), there is a 
larger number of HI than optical galaxies for Eio less than 1 
gal Mpc~ 2 . On running a Kolmogorov-Smirnov (K-S) test, 



3 http:/ /www-wfau. roe. ac.uk/6dfgs/form. html 



vel < 4000 km/s & M H| > 3x1 9 M q| 




Figure 4. Mhi versus Star formation rate for volume limited 
sample vel < 4000 km s _1 and Mhi > 3 x 10 9 Mq for IRAS 
infrared data. Least-squares line of best fit, as described in Fig. 
131 are fitted. Mhi increases with increasing SFR. 

(right panel Fig. |5J, we find the null hypothesis can be re- 
jected with a significance level of 0.1 per cent, therefore the 
two distributions are not the same. At Eio=5 gal Mpc~ 2 the 
difference between the two distributions starts to increase. 
This indicates an increase in the number of optically galaxies 
compared to HI galaxies. We conclude, as Eio increases, the 
number of HI galaxies decreases and the number of optically 
detected galaxies increases. 

4.3 SFR-Ek, Relationship 

We use the two HOPCAT target galaxy lists, as described 
in Sections 13 . 21 and 13 . 1 1 and Table Q to investigate the rela- 
tionship between SFR and Eio for HI galaxies. 

Fig. HJ shows the relationship between SFR and Eio 
for, IR (top panels) and 1.4-GHz (bottom panels) SFRs, for 
galaxies with 1000<vel<6000 km s _1 . Least-squares line of 
best fit, as described earlier, are fitted. The SFR uncertainty 
values are shown as red (min) and black (max). The SFR 
uncertainties for SUMSS are red (min) and black (max) and 
for N VSS, magenta (min) and cyan (max). 

In lGomez et al. f200St and lLewis et al. (2002T) the aver- 
age SFR, for a given density range, decreases with increasing 
Eio. We use scatter plots to compare the SFR and Eio for 
each HI galaxy. In Fig. |S] we do not see a clear decrease in 
the SFR with increasing Eio. No conclusions can be drawn 
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Figure 3. Mjy; versus SFR. Least-squares line of best fit, calculated using the logs of M-hi and SFR. The gradient uncertainty lines are 
3(7. Left panel: IRAS infrared data with the minimum (red) and maximum (black) uncertainties. Right panel: SUMSS & NVSS 1.4-GHz 
radio data with the SUMSS minimum (red) and maximum (black) and NVSS minimum (magenta) and maximum (cyan) uncertainties 
are shown. M# j increases with increasing SFR. 

[b!] 




Figure 5. HOPCAT and 6dFGSDR2 Sio calculated using 6dF2MASS as the background catalogue. Left Panel: Direct comparison of 
the number of HOPCAT and 6dFGSDR2 galaxies with respect to Eio. Right panel: K-S test showing the two distributions are not the 
same. Note at Eio=5 gal Mpc — 2 there is an increase in the difference between the two distributions, indicating an increase in the number 
of optical galaxies. Note the Eio horizontal axis are limited to 40 gal Mpc - 2 



from Fig. |^1 as multiple lines of best fit are possible and the 
flux-limited nature of HOPCAT maybe affecting the results. 

To reduce the flux-limiting effect, we re-analyse SFR 
versus Eio, separated into velocity ranges of 1000 km s . 
In Figs. □ and H (IR) and Figs. M and d (1.4-GHz), the 
SFR per HI galaxy with respect to Eio remains remarkably 
constant. The SFR versus Eio average gradients are; for IR 
0.11±0.05, 0.09±0.15 (Figs.|7|and|HJ and 1.4-GHz 0.12±0.15 
and 0.06±0.15 (Figs. Inland [TO. 

To determine if the relationship is indeed constant, we 



conduct T-tests, where two normally distributed popula- 
tions are equal if the null hypothesis cannot be rejected at 
a low significance level. From the 4 data sets separated into 
5 velocity groups, 18 out of a total of 20 the null hypothesis 
cannot be rejected to a 0.1 per cent significance level. We 
conclude that the SFR for HI galaxies between 1000 and 
6000 km s _1 remains constant regardless of the local galaxy 
density. 
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Figure 6. SFR versus Ei . Top Left panel: HOPCAT-SCosGC-IR. Top Right panel: HOPCAT-6dF2MASS-IR. Bottom Left panel: 
HOPCAT-SCosGC-1.4-GHz. Bottom Right panel: HOPCAT-6dF2MASS-1.4-GHz. The uncertainty values for IR SFRs are shown as 
minimum (red) and maximum (black). The SUMSS uncertainties are red (min) and black (max) and the NVSS uncertainties are 
magenta (min) and cyan (max). Due to the multiple line of best fit possible for these plots no conclusion can be drawn. 



4.4 SFE Versus Em 

We investigate the efficiency of Hi's conversion into stars 
for local galaxy density using SFE and E 10 . Fig. [HI shows 
the relationship between SFE and Eio for IR and Fig. ll2l for 
1.4-GHz SFEs. Again the sample only contains galaxies with 
velocities between 1000 and 6000 km s _1 . Least-squares line 
of best fit, as described earlier, are fitted. The IR ASFE's are 
shown in red (min) and black (max). The 1.4-GHz SFE's un- 
certainty are shown in red (min) and black (max) for SUMSS 
and magenta (min) and cyan (max) for NVSS. For both IR 
and 1.4-GHz, the SFEs remains unchanged with respect to 
Eio- 



To ensure the flux-limited nature of our data is not af- 
fecting the results, we separate our data into velocity ranges 
of 1000 km s _1 . Figs. Hoi and ITU shows SFE versus Eio for 
IR and Figs. [TE] and ED for 1.4-GHz. The line of best fit 
along with its 3<r gradient uncertainty lines are shown. The 
SFE versus Eio average gradients are; IR 0.11T0.05 and 
0.09T0T5 (Figs. E3 and E3 and 1.4-GHz 0T2T0T5 and 
0.06T0T5 (Figs. Upland [T6l. 

T- Tests are conducted to determine if the constant SFE 
for Eio is true. From the 4 data sets separated into 5 veloci- 
ties groups, 19 out of a total of 20 the null hypothesis cannot 
be rejected to a 0.1 per cent significance level. 
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Table 3. T-test results for SFR and SFE versus Si . 



Velocity 


HOPCAT- 


HOPCAT- 


HOPCAT- 


HOPCAT- 


Groupings 


SCosGC-IR 


6dF2MASS-IR 


SCosGS-1.4GHz 


6dF2MASS-1.4GHz 


km s~ x 


P-values 


P-values 


P-values 


P-values 


SFR vs Sio 


Fig.Q 


Figs.lU 


Figs. 


Fig. HOI 


1000<vel<2000 


0.4790t 


0.8710t 


0.3450* 


0.1830* 


2000<vel<3000 


0.0953+ 


0.6740+ 


0.7640* 


0.4410* 


3000<vel<4000 


0.7920+ 


0.8730+ 


0.8940* 


0.4850* 


4000<vel<5000 


0.1500+ 


0.0045+ 


0.9260* 


0.3140* 


5000<vel<6000 


0.3370t 


0.0198t 


0.7750* 


0.4460* 


SFE vs Sio 


Fig. EH 


Figs, rm 


Figs. [TBI 


Fig. [TBI 


1000<vel<2000 


0.2360* 


0.9850 1 " 


0.1990+ 


0.1660* 


2000<vel<3000 


0.4230+ 


0.3470t 


0.4260* 


0.2640* 


3000<vel<4000 


0.8910+ 


0.9700+ 


0.8960* 


0.7940* 


4000<vel<5000 


0.5460+ 


0.3270t 


0.0070* 


0.3420* 


5000<vel<6000 


0.5580t 


0.0137t 


0.5400* 


0.2010* 



*: h=0 - two normally distributed populations are equal if the null hypothesis cannot be rejected at a significance level of 0.1%. 
*: h=l - two normally distributed populations are not equal if the null hypothesis can be rejected at a significance level of 0.1%. 



In Figs. UJl[T3|[Tni and Ell for a single density, the SFE 
varies, however as Eio increases the SFE remains unchanged. 
Our conclusion is that the efficiency, with which the HI 
present in galaxies is converted to stars, remains constant 
regardless of local surface density. 



5 DISCUSSION 

In this section we discuss the relationships between Mhi 
and SFR, SFR and Eio, and SFE and Eio. We also compare 
optical and HI detected galaxy samples. From these results 
we determine which of our two theories explain the SFR- 
density relation. 

5.1 HI mass-SFR Relationship 

We investigate the relationship between the SFR and Mhi 
HOPCAT galaxies with IR and 1.4-GHz SFRs. We find 
Mhi increases with increasing SFR. This result is in general 
agreement with the Schmidt law, a relationship between a 
galaxy's gas surface density and the disk average SFR as 
disc ussed inlKennicutt fl99&t) . 

iMirabel fc Sanders (198a) . investigating 92 luminous 
infrared galaxies find that 88 galaxies have M^j in the range 
of 5 x 10 8 to 3 x 1O 1O M . They conclude the FIR luminosi- 
ties have only a loose correlation with the Mhi- However, 
our study looks at HI galaxies and calculates SFR from IR 
and 1.4-GHz fluxe s and finds a s trong correlation between 
the SFR and Mh r. lStevens f20 05'l also investigated the rela- 
tionship between Mhi and SFR using a volume limited sub- 
sample of HICAT, containi ng 85 FIR and 277 radio (NVSS 
and SUMSS) based SFRs. IStevens (20051) also finds M H i 
increases with increasing SFR. 

Although we do not know where star formation occurs 
in a galaxy with respect to the HI detection, the correlation 
between the Mhi and SFR shows HI is an indication of star 
formation. In fact, as our SFRs are from IR and 1.4-GHz 
fluxes, which are indicators of current star formation, the 



Mhi and SFR correlation suggests HI indicates 'current' 
star formation in galaxies. 

5.2 HI radio and Optical galaxy samples 

Our result from investigating the difference between opti- 
cal and HI detected galaxy samples show there is a much 
larger population of elliptical galaxies than HI detected 
galaxies with increasing local surface density. This result 
agrees with stu dies into the morphology-density relation 
iDressler f 1 98011 etc.), which find for an increasing local 
galaxy density there is a decrease in the number of spiral 
galaxies and an increase in the number of elliptical galaxies. 
It is expected that this decrease in the number of current 
star forming HI galaxies with increasing Eio will have a di- 
rect effect on the results from any SFRs versus Eio analysis 
for optical samples. 

5.3 SFR and E i0 

From our investigation into SFR and Eio we conclude the 
SFRs for our HI galaxies are not affected by increasing Eio. 

As HOPCAT is based on the Parkes HI All Sky Sur- 
vey (HIPASS) there may be some concern about the num- 
ber of low surface brightnes (LSB) galaxies in the sample. 
The HIPASS galaxy sample is dominated by late-type galax- 
ies as is expecte d due to the HI selecti on (e.g. Helmboldt 
et al. 2004). As iHelmboldt et al. (20041 show, Hl-selected 
galaxies tend to have lower luminosity and surface bright- 
ness than a comparison sample of optically-selected galaxies. 
This difference is not so great for our HOPCAT sample be- 
cause our mean redshift is much greater than the Helmboldt 
sample. Our median luminosity is Mb = —18.8 mag com- 
pared to Mb = —18.5 mag in their Hl-selected sample and 
Mb = —19.1 mag in their optical sample. Some 15 per cent 
of their Hl-selected galaxies have LSB compared to the op- 
tical sample: the fraction of LSB galaxies in HOPCAT will 
be smaller because of the higher median luminosity. If we 
consider the final HOPCAT sample with reliable densities 
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100 




Figure 7. SFR versus Sio for HOPCAT-SCosGC-IR, separated into velocity ranges of 1000 km s _1 . Least -squares line of best fit, 
calculated using the logs of Sio and SFR. The gradient uncertainty lines are 3a. The SFR remains unchanged with increasing £io- We 
compare the distributions of low and high Sio using T-tests and find the null hypothesis cannot be rejected to a 0.1 per cent significance 
level. See Table [3] for T-test results. 



and star formation rates used in our analysis, only 1.3 per 
cent of the galaxies are fainter than the Mb, =-17 mag faint 
absolute magnitude limit applied in the density calculations. 

In lGomez et al. f 20031) and lLewis et al. f 2002ft the de- 
crease in the average SFR with an increase in local sur- 
face density, for a given density range, is explained by 
their optical sample. An optical sample contains both spi- 
ral and elliptical galaxies. For areas of increasing local 
galaxy density, such as compact galaxy groups and poor 
and rich clusters, the number of spiral galaxies decreases 
as the number of elliptical gala xies increases. S tudies into 
the morphology-density relation jDressler etc.), show 



this change in morphology type for increasing local galaxy 
density. Elliptical galaxies have lit tle or no current star for- 
mation, hence the average SFR, as lGomez et al. f2003ft and 
iLewis et, al. (20021) use, will decrease in areas of increasing 
local surface density where elliptical galaxy numbers in- 
crease and spiral galaxy numbers decrease. As our sample 
does not have elliptical galaxies and the SFRs are for galax- 
ies with current star formation, a more direct picture of SFR 
with respect to local surface density is seen. We clearly show 
the HI galaxies' SFRs are relatively constant regardless of 
local surface density. 

iGomez et al. f2003ft and lLewis et al.' f2002ft describe a 
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Figure 8. SFR versus Eio for HOPCAT-6dF2MASS-IR, separated into velocity ranges of 1000 km s -1 . Least -squares line of best fit, 
as described in Fig. |7| are fitted. Again the SFR remains unchanged with increasing Eio- We compare the distributions of low and high 
Sio using T-tests and find the null hypothesis cannot be rejected to a 0.1 per cent significance level. See Table 131 for T-test results. 



characteristic Eio break at 1 gal Mpc -2 where the average 
SFR decreases sharply. As our scatter plots are a direct com- 
parison of SFR per HI galaxy for each galaxies local surface 
density, using lines of best fit, no break is seen. However, at 
5 gal Mpc -2 (right panel of Fig.|KJ, the increasing gap indi- 
cates a decrea se in the number of HI g alaxies relative to op - 
tical galaxies. iLewis et al. (2002ft and|G omcz et al. (2003) 
use a faint absolute magnitude limit of -19 mag and we use 
-17 mag which results in our larger Eio values. Our decrease 
in H I galaxy numbers a t Ei o=5 gal Mpc~ 2 m ay coincide 
with lGomez et al. (20031 and lLewis et al. ( 2002) character- 
istic Sio break at 1 gal Mpc -2 . 

IStevens (20051) . also using HICAT, finds the SFR for 
galaxies in groups is smaller than field galaxies. This differs 



from the results presente d in this paper for the same sample. 
However, IStevens (20051) uses a group finder algorithm to 
find HI galaxies in groups (143) and compares their SFR to 
those for field g alaxies (17). A s the resulting density is 3D 
and the analysis [Stevens (2005ft uses is significantly different 
from ours, no direct comparis on is possible. 

However, one observation IStevens (2005) makes is there 
is a larger population of red galaxies in groups than in the 
field. He concludes this could partially explain the lower SFR 
in group galaxies, which agree with our results. It is clear 
the lo wer average SFR in regio ns of high local surf ace den- 
sity in ILewis et al. (2002ft and iGomez et al. (2003) can be 
attri buted to their large r population of elliptical galaxies. 

iBaloeh ' et al. (2004) using Ha equivalent widths to in- 
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Figure 9. SFR versus Sio for HOPCAT-SCosGC-1.4GHz, separated into velocity ranges of 1000 km s _1 . Least -squares line of best fit, 
as described in Fig. [7| are fitted. The SFR remains unchanged with increasing Sio. We compare the distributions of low and high Sio 
using T-tests and find the null hypothesis cannot be rejected to a 0.1 per cent significance level. See Table 131 for T-test results. 



vestigate the population of star-forming galaxies in SDSS 
and 2dFGRS and finds 'the distribution of line strength' for 
these star- forming galaxies is independent of environment. 
This appears to be in agreement with our results. 

Our conclusion still remains, for HI galaxies between 
1000 and 6000 km s~ x in the southern sky, the SFR does 
not depend on the local galaxy density. 



5.4 Fewer Star Forming Galaxies or Suppressed 
Star Formation? 

In Section f5. II we show a strong correlation exists between 
Mhi and SFR leading to the conclusion our HI galaxies dis- 



play current star formation. As discussed in Section r5.2l the 
number of HI galaxies decreases sharply in regions of high 
local surface densit y. Th e difference between this work and 
iGomez et al. ( 20051) and iLewis et al. f2002) is our sample 
contains current star forming galaxies and theirs are optical 
galaxy samples which contain galaxies with little or no cur- 
rent star formation. We also find in Section 15.31 the SFRs 
do not depend on Em. Our results in Section [4.41 show the 
efficiency, with which HI is converted to stars, is not affected 
by an increase in the local surface density. We now conclude 
which theory explains the SFR-density relation. 

We find there are, indeed, fewer star forming (HI) galax- 
ies in regions with high local surface density, which affects 
the average SFR for optical galaxy samples. As the SFR and 



12 M. T. Doyle et al. 




Figure 10. SFR versus Eio HOPCAT-6dF2MASS-1.4GHz, separated into velocity ranges of 1000 km s _1 . Least-squares line of best fit, 
as described in Fig. [7| are fitted. The SFR remains unchanged with increasing Eio. We compare the distributions of low and high Eio 
using t-tests and find the null hypothesis cannot be rejected to a 0.1 per cent significance level. See Table [3] for T-test results. 



SFE do not depend on the local surface density, the SFR- 
density relation is due to the low number of star forming 
(HI) galaxies, in regions of high local surface density and is 
not due to the suppression of star formation in these regions. 



6 SUMMARY 

We use an optically matched HI detected galaxy catalogue, 
HOPCAT, to investigate SFR, SFE and the local surface 
density to explain the SFR-density relation. In regions of 
high local galaxy density either, there are fewer star form- 
ing galaxies, or the galaxies capable of forming stars are 
present but some physical process is suppressing their star 



formation. The local surface density, Eio, is calculated us- 
ing HOPCAT with SCosGS and 6dF2MASS as background 
catalogues. IRAS, SUMSS and NVSS fluxes are used to 
calculate the SFRs and SFEs. For the SFR-M/// inves- 
tigations we use 2 different data sets using 6dF2MASS 
with IR and 1.4-GHz SFRs, HOPCAT-6dF2MASS-IR (1405 
galaxies) and HOPCAT-6dF2MASS-1.4GHz (561 galax- 
ies). For the SFR-SFE-Eio investigations we have 4 dif- 
ferent data sets using SCosGS and 6dF2MASS each with 
IR and 1.4-GHz SFRs, HOPCAT-SCosGS-IR (923 galax- 
ies), HOPCAT-6dF2MASS-IR (1180 galaxies), HOPCAT- 
SCosGS-1.4GHz (339 galaxies) and HOPCAT-6dF2MASS- 
1.4GHz (478 galaxies). We find: 

1. A strong correlation exists between Mhi and SFR 
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Figure 11. SFE versus Ei . Left panel: HOPCAT-SCosGC-IR. Right panel: HOPCAT-6dF2MASS-IR. Least-squares line of best fit, 
calculated using the logs of Sio and SFE, with 3<r gradient uncertainty lines. The minimum (red) and maximum (black) uncertainties 
are shown. The SFE remains unchanged with increasing local galaxy density. 




Figure 12. SFE versus Ei . Left panel: HOPCAT-SCosGC-1.4-GHz. Right panel: HOPCAT-6dF2MASS-1.4-GHz. Least-squares line of 
best fit, as described in Fig. 1111 are fitted. The SFE remains unchanged with increasing Sio. The SFE's uncertainty are shown in red 
(min) and black (max) for SUMSS and magenta (min) and cyan (max) for NVSS. Like the infrared, the SFE is unchanged with respect 
to local galaxy density. 



for galaxies out to 10000 km s _1 in the southern sky. The 
SFRs are calculated using IR and 1.4-GHz fluxes, both of 
which are indicators of current star formation. We conclude 
HI galaxies display current star formation and their Mhi 
and SFR are strongly related. 

2. The number of HI galaxies decrease s with increasin g 
local surface density. This agrees with iDressler ( 198(3) . 



[ Hashimoto et al. (lflflish . iMartmez et al. C200S ), 

[ Goto et al. (20031). ISodre fc Mateus (200; i. 

iNuiiten et al. (2005T) and IChristlein fc Zabludoff (20051) . 
where optical samples show an increase in the number of 
non-current star forming elliptical galaxies with a decrease 
in current star forming galaxy nu mbers for increasing local 
galaxy density. This is also seen in Balo gh et al. (2 004). 
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Figure 13. SFE versus Sio for HOPCAT-SCosGC-IR, separated into velocity ranges of 1000 km s _1 . Least -squares line of best fit, as 
described in Fig, 1111 are fitted. The SFE remains unchanged with increasing £10. We compare the distributions of low and high Sio 
using T-tests and find the null hypothesis cannot be rejected to a 0.1 per cent significance level. See Table 131 for T-test results. 



3. For our HI galaxies, the SFR does not depend on lo- 
cal surface density. Any decrease in the average SFR with 
increasing local surface density for optical galaxy samples, 
is explained by the decrease in the number of current star 
forming galaxies and an increase in elliptical galaxies in re- 
gions with increasing local surface density. 

4. The efficiency with which the HI present in the HI 
galaxies is converted into stars is the same regardless of local 
surface density. 

5. The SFR-density relation is due to a decrease in the 
number of star forming (HI) galaxies in regions of high local 
(surface) galaxy density and not due to the suppression of 
star formation in these regions. 
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Figure 14. SFE versus Sio for HOPCAT-6dF2MASS-IR, separated into velocity ranges of 1000 km s _1 . Least -squares line of best fit, 
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